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ABSTRACT 

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z < 0.5) detected at high signal-to- 
noise in the first Planck all-sky data set. The sample spans approximately a decade in total mass, 2 x lO''' Mq < M500 < 2 x 10" Mq, where M500 
is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray 
data, we investigate the relations between Ysqq, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass Mj ^oo, 
temperature Tx, luminosity Lx,500> SZ signal analogue Yx^sm = 500 x Tx, and total mass M^m- After correction for the effect of selection bias on 
the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived 
from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for 
all-sky statistical studies of galaxy clusters. 

Key words. Cosmology: observations. Galaxies: cluster: general. Galaxies: clusters: intracluster medium. Cosmic background radiation. X-rays: 
galaxies: clusters, Planck satellite 

1. Introduction (o-j/meC^) J Pdl. Here P oc r is the ICM thermal electron 

pressure, where is the density and T is the temperature, crj 

The X-ray emittmg gas m galaxy clusters mduces mverse Thomson cross section, is the electron rest mass, and 

Compton scattermg of Cosmic Microwave Background (CMB) ^ ^j^^ ^p^^^ jigj^^ S2 signal integrated over the clus- 

photons, shifting their frequency distribution towai-ds higher en- ^„ proportional to the integrated Compton parameter, 

ergies. First discussed in 1972 by Sunyaev & Zeldovich, the ^^^^ ^^^^ ^2 ^ (crr/m^c^) [PdV, where Da is the angular 

scattering produces a characteristic distortion of the CMB spec- j- » » ii, 

distance to the source 

trum in the direction of a cluster known as the thermal Sunyaev- 
Zeldovich (SZ) effect. It is directly proportional to the Compton Clusters are currently thought to form via the hierarchical 
parameter y, a measure of the thermal electron pressure of the gravitational collapse of dark matter haloes, so that their num- 
intracluster medium (ICM) gas along the line of sight: y - ber as a function of mass and redshift is a sensitive indicator 

of the underlying cosmology. The ICM is formed when gas falls 

* Corresponding author: G.W. Pratt, gabriel . prattOcea . fr into the dark matter gravitational potential and is heated to X-ray 
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emitting temperatures by shocks and compression. The scale- 
free nature of this process implies that simple power law rela- 
tionships exist between the total halo mass and various other 
physical properties (e.g., Bertschinger 1985; Kaiser 1986) such 
as X-ray temperature T or luminosity L (e.g., Voit 2005; Arnaud 
et al. 2005, 2007; Pratt et al. 2009; Vikhlinin et al. 2009). As 
the total mass is not directly observable, such mass proxies are 
needed to leverage the statistical power of various large-scale 
surveys for cosmological applications. Since the gas pressure is 
directly related to the depth of the gravitational potential, the 
quantity Fsz is expected to scale particularly closely with the 
total mass, a claim supported by recent numerical simulations 
(e.g.. White et al. 2002; da Silva et al. 2004; Motl et al. 2005; 
Nagai 2006; Wik et al. 2008; Aghanim et al. 2009). SZ surveys 
for galaxy clusters thus have great potential to produce compet- 
itive cosmological constraints. 

In a few short years, SZ observations have progressed from 
the first spatially resolved observations of individual objects 
(Pointecouteau et al. 1999, 2001; Komatsu et al. 1999, 2001), to 
the first discoveries of new objects (Staniszewski et al. 2009), to 
large-scale survey projects for cosmology such as the Atacama 
Cosmology Telescope (Kosowsky 2003, ACT) and the South 
Pole Telescope (Carlstrom et al. 2011, SPT). Indeed, first cos- 
mological results from these surveys have started appearing 
(Vanderlinde et al. 2010; Sehgal et al. 2011). Attention is now 
focussing on the shape and normalisation of the pressure profile 
(e.g., Nagai et al. 2007; Arnaud et al. 2010; Komatsu et al. 201 1 ), 
calibration of the relationship between D\ Fsz and the total 
mass for cosmological applications (e.g., Bonamente et al. 2008; 
Marrone et al. 2009; Arnaud et al. 2010; Melin et al. 2011), 
comparison of the measured SZ signal to X-ray predictions 
(Lieu et al. 2006; Bielby & Shanks 2007; Afshordi et al. 2007; 
Komatsu et al. 2011; Melin et al. 2011), and the relationship 
between Isz and its X-ray analogue Fx.soo (e.g., Andersson 
et al. 2010). First introduced by Kravtsov et al. (2006), the latter 
quantity is defined as the product of Mg, the gas mass, and Tx, 
the spectroscopic temperature excluding the core regions. As the 
link between Fx.soo and Y^z depends on the relationship be- 
tween the gas mass weighted and X-ray spectroscopic tempera- 
tures, it is a sensitive probe of cluster astrophysics. 

In the following, we use a subsample of Planck^ Early 
Release Compact Source Catalogue SZ (ESZ) clusters, consist- 
ing of high signal-to-noise ratio Planck detections with deep 
XMM-Newton archive observations, to investigate the local (z < 
0.5) SZ scaling relations. Given its all-sky coverage and high 
sensitivity, Planck is uniquely suited to this task, allowing high 
signal-to-noise ratio detection of many hot, massive systems 
that do not appear in other SZ surveys due simply to their lim- 
ited area; correspondingly, the large field of view and collecting 
power of XMM-Newton make it the ideal instrument to observe 
these objects in X-rays out to a significant fraction of the virial 
radius. Here we investigate the relationship between SZ quan- 
tities and X-ray quantities, making full use of the exceptional 
quality of both data sets. Two complementary companion papers 
(Planck Collaboration 201 If and Planck Collaboration 201 Ih) 
harness the statistical power of the Planck survey by analysing 
the SZ flux-X-ray luminosity and SZ flux-optical richness re- 



' Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



lations, respectively, using a bin-averaging approach. Two fur- 
ther companion papers present the parent catalogue (Planck 
Collaboration 201 Id) and XMM-Newton validation observations 
of newly-discovered clusters (Planck Collaboration 201 le). 

In this paper we adopt a ACDM cosmology with Hq = 
70 km s ' Mpc-i, Qm = 0.3 and Qa = 0.7. The factor E(z) = 
is the ratio of the Hubble constant at red- 
shift z to its present day value. The variables M500 and /^soo are 
the total mass and radius corresponding to a total density con- 
trast 500pc(z), where pdz) is the critical density of the Universe 
at the cluster redshift; thus M500 = (47r/3)500pc(z)/?5o,). The 
quantity ix.soo is defined as the product of Mgsoo, the gas 
mass within Rsoq, and Tx, the spectroscopic temperature mea- 
sured in the [0.15-0.75] R5Q0 aperture. The SZ signal is de- 
noted Fsoo throughout. This quantity is defined by F500 = 
(a-j/nieC^) J PdV. Here Da is the angular distance to the sys- 
tem, cTj is the Thomson cross-section, c the speed of light, 
the electron rest mass, P ocn^T is, the pressure (the product of the 
electron number density and temperature), and the integration is 
performed over a sphere of radius Rsm- The quantity Ksoo is 
the spherically integrated Compton parameter, and Y^oq is pro- 
portional to the flux of the SZ signal within /^soo. 

2. The ESZ catalogue and the Planck-XMM-Newwn 
archive subsample 

2.1. Planck and the ESZ Catalogue 

Planck (Tauber et al. 2010; Planck Collaboration 2011a) is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31' to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; Bersanefli et al. 2010; 
Mennella et al. 2011 covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolometers 
cooled to 0. 1 K. Polarization is measured in all but the highest 
two bands (Leahy et al. 2010; Rosset et al. 2010). A combina- 
tion of radiative cooling and three mechanical coolers produces 
the temperatures needed for the detectors and optics (Planck 
Collaboration 2011b). Two Data Processing Centres (DPCs) 
check and calibrate the data and make maps of the sky (Planck 
HFI Core Team 2011b; Zacchei et al. 2011). Planck's sensitiv- 
ity, angular resolution, and frequency coverage make it a pow- 
erful instrument for galactic and extragalactic astrophysics as 
well as cosmology. Early astrophysics results are given in Planck 
Collaboration (201 le-u) . 

The basic data set used in the present paper is the Planck 
Early Release Compact Source Catalogue SZ (ESZ) sample, de- 
scribed in detail in Planck Collaboration (201 Id). The sample 
is derived from the highest signal-to-noise ratio detections (S/N 
> 6) in a blind multi-frequency search in the all-sky maps from 
observations obtained in the first ten months of the Planck mis- 
sion. 

2.2. The Planck-XMM-Newton archive subsample 

Cross-correlation of the ESZ subsample with the Meta 
Catalogue of X-ray Clusters (MCXC; Piffaretti et al. 2010) pro- 
duced 158 matches with known X-ray clusters. As shown in 
Fig. 1, these objects lie at a redshift z < 0.5 and cover approx- 
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Fig. 1: The 158 Planck ESZ clusters already observed in X-rays. 
Masses are estimated from their X-ray luminosity as described in 
PifFaretti et al. (2010). The 62 clusters analysed and discussed in this 
paper are indicated by the red squares. The dashed line represents the 
locus at which i?5oo ~ 12'. 



imately a decade in mass-. A search for these clusters in the 
XMM-Newton Science Archive ' produced a combined Planck- 
XMM-Newton archive sample of 88 objects as of July 2010, in- 
dicated by blue crosses in Fig. 1 . 

As detailed below in Sect. 3, we used different X-ray data 
processing techniques depending on cluster angular extent. More 
specifically, if the source extent lies well within the XMM- 
Newton field of view then the X-ray background can be charac- 
terised using a source-free region of the observation, while clus- 
ters with a larger angular extent require simultaneous source and 
background modelling. We label these classes of clusters as A 
and B, respectively. The dashed line in Fig. 1 illustrates the ra- 
dius at which /?5oo < 12', corresponding to the maximum angular 
extent within which the X-ray background can be characterised 
in a single XMM-Newton field of view. Using this criterion, we 
divide the Planck-XMM-Newton archive sample into 58 A clus- 
ters and 30 B clusters. 

Not all of the clusters in the full Planck-XMM-Newton 
archive sample are used in the present paper. Some observa- 
tions in the A cluster list were excluded because soft proton 
solar flare contamination had rendered the observations unus- 
able, or because the object had not yet been observed at the time 
of the archive search, or because the target was a clear multi- 
ple system unsuited to a spherically-symmetric analysis. For the 
B clusters, in addition to the high-luminosity systems already 
published in Bourdin & Mazzotta (2008), we prioritised those 
where the XMM-Newton field of view was expected to cover the 
largest possible fraction of /?5oo, corresponding to objects with 
the lowest estimated mass in Fig. 1. The final sample of 62 sys- 
tems consists of 44 A objects and 18 B objects. While the sample 
is neither representative nor complete, it represents the largest, 
highest-quality SZ-X-ray data set currently-available. 



^ Estimated from the X-ray luminosity-mass relation of Pratt et al. 
(2009), as detailed in Pitfaretti et al. (2010). 
^ http://xmin.esac.esa.int/xsa/ 



3. X-ray cluster properties 

3. 1. X-ray data processing: A clusters 

For the A clusters, we use the latest version (vlO.O) of the XMM- 
Newton-SAS, ensuring that the most recent calibration correc- 
tions are applied to the X-ray photons. Event lists are processed, 
cleaned for periods of high background, PAXXERN-selected and 
corrected for vignetting as detailed in Pratt et al. (2007). 

The quiescent XMM-Newton background is dominated by 
events due to charged particles. We subtract this component 
by using a background file built from stacked observations ob- 
tained with the filter wheel in the CLOSED position, recast to the 
source position and renormalised using the count rate in the high 
energy band free of cluster emission"^. The remaining compo- 
nents are due to the cosmic X-ray background produced by unre- 
solved sources and a diffuse soft X-ray contribution attributable 
to the Galaxy. The A sample is selected so that R^qo < 12', al- 
lowing us to model these remaining components using emission 
from an annular region external to the cluster emission as de- 
tailed in Croston et al. (2008) and Pratt et al. (2010). 

Point sources were identified from the small scales of 
wavelet-decomposed images in the [0.3-2] and [2-5] keV bands. 
These sources were excluded, with the exclusion radius matched 
to the variation of the PSF size across the detector. We also 
masked any well-defined substructures that produced prominent 
secondary maxima and were visible in the larger scales of the 
wavelet decomposition process. 

Surface brightness profiles were extracted from the data in 
3'.'3 bins centred on the X-ray peak. Finally, a non-parametric 
regularisation method was used to derive the deprojected, PSF- 
corrected density profiles, «e(r), as described in Croston et al. 
(2008). 



3.2. X-ray data processing: B clusters 

For each object in the B cluster sample, a merged energy- 
position photon cube was built from the various observations 
of a given target. The cube was built from soft proton-cleaned 
events from each camera, generated with vlO.O of the XMM- 
Newton-SAS, to which an effective exposure and a background 
noise array were added. The exposure array was computed from 
the effective exposure time, with corrections for spatially vari- 
able mirror effective areas, filter transmissions, CCD pixel area, 
chip gaps and bad pixels, as appropriate. The background noise 
array was modelled as the sum of components accounting for 
the Galactic foreground and cosmic X-ray background, plus 
charged particle-induced and out-of-time events. Full details of 
the method are given in Bourdin & Mazzotta (2008). 

The Galactic foreground is a critical model component in 
the case of the B clusters. These objects often extend over the 
full XMM-Newton field of view, so that the cluster emission can- 
not be spatially separated from the foreground components. We 
thus constrained the foreground components using a joint fit of 
cluster emissivity and temperature in an external annulus cor- 
responding to ~ RsQQ- Despite the degeneracy of this estimate 
with the cluster emissivity itself, in all cases the temperature ob- 
tained in this annulus was found to be lower than the average 
cluster temperature, as is commonly observed in clusters allow- 
ing full cluster-foreground spatial separation (e.g., Pratt et al. 
2007; Leccardi & Molendi 2008). 



We excluded a 5' region around the cluster centre to avoid contam- 
ination from residual cluster emission. 



3 



Planck Collaboration: Local galaxy cluster SZ scaling relations 




10^* 10^^ 0.01 0.10 1.00 

M500 [Md Radius [R500] 

Fig. 2: Left: Lx,5oo - M500 relation of the Planck-XMM-Newton archive sample compared to REXCESS, a representative X-ray cluster sample. 
Luminosity is estimated interior to i?5oo in the [0.1-2.4] keV band, and mass from the M500 - ix.soo relation of Arnaud et al. (2010). The solid 
red line is the fit to the REXCESS sample only (Pratt et al. 2009). Right: Scaled density profiles of the 62 systems in the Planck-XMM-Newton 
archive sample. Profiles have been corrected for projection and PSF effects as described in the text. Systems classified as cool cores are indicated 
by blue lines. 



The ICM density profiles in the B cluster sample were then 
derived using the analytic distributions of ICM density and tem- 
perature introduced by Vikhlinin et al. (2006). These parametric 
distributions were projected along the line of sight, convolved 
with the XMM-Newton PSF, and fitted to the observed projected 
cluster brightness and temperature profiles. The resulting den- 
sity profiles, n^ir), were used to derive X-ray quantities for each 
cluster as described below. 



3.3. X-ray quantities 

For the current generation of high-resolution X-ray telescopes, 
the effective limiting radius for high-quality nearby observations 
of the type discussed here is Rsoo. Beyond this radius, the effect 
of the variable background becomes dominant and the uncertain- 
ties begin to become difficult to quantify. In addition, as shown 
by Evrard et al. (1996), Rs{H) is also the radius within which clus- 
ters are relatively relaxed. We estimate the X-ray quantities for 
each cluster self-consistently within /?5oo using the M500 - ix,500 
relation given in Arnaud et al. (2010, see also Pratt et al. 2010), 
assuming standard evolution, viz.. 



E(z)-"Msi 



10 



14.567±0.010 



2x lO'^MokeV 



(1) 



The radius R^qo was calculated iteratively as described in 
Kravtsov et al. (2006). Using Eq. 1 and the definition of M500 
and Fx,500, an equation of the form TJ^^g - C[Mg 500 Tx]" must 
be solved. Starting from an initial temperature measurement, the 
equation is solved for Rsa), with the gas mass computed from 
the density profiles discussed above in Sects. 3.1 and 3.2. A 
new temperature is then estimated within [0.15 - 0.75] R5Q0 and 
the procedure is repeated until convergence. Once converged. 



we measure the temperature within an aperture corresponding to 
Rsm and calculate the [0.1 - 2.4] keV and [0.5 - 2] keV band lu- 
minosities as described in Pratt et al. (2009). The resulting X-ray 
quantities are listed in Table 1 . 

The left-hand panel of Fig. 2 shows the Lx,500 - A/500 relation 
of the Planck-XMM-Newton archive sample, where Lx.soo is the 
X-ray luminosity estimated interior to R^m in the [0.1-2.4] keV 
band^, and the mass is estimated from the Mgoo - ix,500 relation 
given in Equation 1 . The data are compared to the equivalent re- 
lation from REXCESS, a sample designed to be representative 
of the X-ray cluster population (Bohringer et al. 2007). One can 
see that the Planck-XMM-Newton archive clusters are all mas- 
sive, luminous systems, as expected for objects detected in SZ 
at high S/N by Planck. They follow the general trend exhibited 
by REXCESS (Pratt et al. 2009), but extend to higher mass and 
luminosity. 

3.4. Scaled gas density profiles and cool core subsample 

The scaled gas density profiles of the full sample of 62 clusters 
are shown in the right-hand panel of Fig. 2, where each profile 
has been corrected for evolution and scaled to R500. As has been 
seen in other cluster samples (e.g., Croston et al. 2008), there is 
a large amount of scatter in the central regions, extending out to 
~ 0. 15 /?5oo, beyond which the profiles rapidly converge. 



' The standard [0. 1 - 2.4] keV band self-similar evolution factor is 
E(zy^. However, the observed scaling relations are steeper than self- 
similar (see e.g., Pratt et al. 2009), leading to a dependence of the evo- 
lution factor on the quantity to which it is applied. In the present work 
we use an evolution factor of E(zy''^^ , appropriate for bolometric lumi- 
nosities, for comparison to REXCESS (see also Maughan 2007). 
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Table 1 : X-ray and SZ properties. 



Name RA Dec z R500 Tx M^joo J'xsoo ^500 Msoo Lx.soo CC 





[deg] 


[deg] 




[kpc] 


[keV] 


[10'" Mo] 


[10'" Mo keV] 


[lO-^Mpc"] 


[10'" Mo] 


[10"" ergs-'] 






RXC J0014.3-3022 


3.58 


-30.38 


0.307 


1358 


7.72± 0.25 


1.65± 0.01 


12.73± 0.51 


1.74± 0.21 


9.78± 0.21 


13.35+ 0.09 




A85 


10.44 


-9.37 


0.052 


1206 


5.78± 0.22 


0.66± 0.01 


3.84± 0.19 


0.47± 0.05 


5.30± 0.31 


4.65± 0.02 


/ 


RXCJ0043.4-2037 


10.84 


-20.61 


0.292 


1152 


5.82± 0.20 


0.88± 0.01 


5.10± 0.20 


1.40± 0.17 


5.88± 0.14 


8.26± 0.08 




A119 


14.02 


-1.30 


0.044 


1114 


5.40± 0.23 


0.45± 0.01 


2.45± 0.14 


0.27± 0.03 


4.12± 0.23 


1.52± 0.01 




RXC J0232.2-4420 


38.06 


-44.37 


0.284 


1223 


6.41± 0.20 


1.07± 0.01 


6.86± 0.26 


0.86± 0.13 


6.95± 0.15 


12.53± 0.09 


/ 


A40I 


44.73 


13.56 


0.075 


1355 


7.26± 0.44 


1.02± 0.04 


7.43± 0.58 


0.83± 0.08 


7.65± 0.67 


5.82± 0.04 




RXC J0303. 8-7752 


46.00 


-77.88 


0.274 


1251 


7.88± 0.36 
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Notes. The temperature Tx is measured in the [0.15 - 0.75] i?5oo region, and the luminosity Lx.500 is measured interior to i?5oo in the [0.1 - 2.4] 
keV band. The mass Mjoo is estimated from the M500 - J'x.soo relation given in Eqn. 1 . The final column indicates whether the cluster is classified 
as a cool core system, defined as described in Sect. 3.4. 



It is well-known that some clusters exhibit so-called cool 
cores, central regions of very dense gas where the cooling time is 
less than the Hubble time (e.g., Jones & Forman 1984). Such ob- 
jects have very high X-ray luminosities and extremely low cen- 



tral entropies that tend to set them apart from the rest of the X-ray 
cluster population (e.g., Fabian et al. 1994; Pratt et al. 2010). In 
addition, the cuiTent consensus is that these systems represent a 
generally more relaxed subset of the cluster population (although 
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see Burns et al. 2008 for a dissenting view). Following Pratt et al. 
(2009), we estimated the central gas density We.o using a /? model 
fit to the gas density profile interior to 0.05 Rsm, and classified 
objects with E{zy^ We.o > 4 x 10"^ cm"-' as cool core systems. In 
total, 22/62 clusters in the present sample are classified as such. 
These are plotted in blue in Fig. 2 and in all following plots. 



4. SZ cluster properties 

4.1. Optimisation of the SZ flux extraction 

The basic SZ signal extraction procedure is described in full 
in Planck Collaboration (2011d). In brief, this procedure con- 
sists of applying multi-frequency matched filters (MMF, Melin 
et al. 2006), that incorporate prior knowledge of the signal, to 
the Planck maps. Specifically, the ICM pressure is assumed to 
follow the universal profile shape derived from the REXCESS 
sample by Arnaud et al. (2010). The SZ flux is computed by in- 
tegrating along the line-of-sight and normalising the universal 
pressure profile. Each profile is truncated at 5 xR^qq , effectively 
giving a measure of the flux within a cylinder of aperture radius 
5 X /?5oo, and then converted to the value in a sphere of radius 
7?50() for direct comparison with the X-ray prediction. This is the 
fundamental SZ quantity used in the present pape/\ and we refer 
to it throughout as Y^oo- 

Section 6 of Planck Collaboration (20 lid) shows that the 
cluster flux derived from blind application of the MMF algo- 
rithm is systematically larger than X-ray expectations. This dis- 
crepancy is a result of overestimation of the cluster size 6500 due 
to the freedom to optimise significance with position and size. 
As shown in Planck Collaboration (201 Id), if the SZ signal is 
instead extracted from a region centred on the X-ray position 
with size ^500 estimated from the X-ray luminosity-mass rela- 
tion, the SZ flux is in better agreement with X-ray expectations. 
When additional constraints on the cluster size are available, the 
SZ flux extraction can be further optimised. 

With the present cluster sample we can make use of the 
higher-quality estimate of the X-ray size 6500, derived from /?5oo, 
measured using the Msoo - Fx.soo relation as detailed in Sect. 3.3. 
Appendix A details the improvement in SZ flux extraction when 
these higher-quality size estimates are used. For each cluster in 
the sample, we thus re-ran the SZ flux extraction, calculating 
Fsoo with the X-ray position and size fixed to the refined values 
derived from the high-quality XMM-Newton observation. 



4.2. Robustness tests specific to iocai sample 

Section 6 of Planck Collaboration (201 Id) details various ro- 
bustness tests relevant to all Planck SZ papers, including inves- 
tigation of the cluster size-flux degeneracy discussed above, the 
impact of the assumed pressure profile used for cluster detec- 
tion, beam-shape effects, colour corrections, contamination by 
point sources, and discussion of the overall error budget. For 
the present sample we undertake two further robustness tests: 
the first is related to the impact of radio source contamination; 
the second examines the impact of the assumed pressure profile 
shape on the derived Y^qq. 
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Fig. 3: Comparison of Ysoo from extraction using the baseline pressure 
profile with that from the best fitting REXCESS cool core and morpho- 
logically disturbed pressure profiles given in Appendix C of Arnaud 
et al. (2010). Cool core systems are marked as blue stars, other sys- 
tems as black dots. Top: universal pressure profile vs cool-core pressure 
profile; Bottom: universal pressure profile vs morphologically-disturbed 
pressure profile. The trend with Rsqq is due to the inability of the Planck 
beam to resolve different profiles at small angular size. The effect is 
small (maximum ~ 10 per cent) and quasi-symmetric (as expected), so 
no bias is introduced. 



4.2.1 . Contamination by point sources 

Contamination by point sources can affect extraction of the SZ 
parameters, and have implications for astrophysical studies of 
clusters or further cosmological applications (Aghanim et al. 
2005; Douspis et al. 2006). We have thus checked the possible 
effect of radio galaxies on the derived Ksoo by combining data 
from SUMSS (Bock et al. 1999, a catalogue of radio sources 
at 0.85 GHz), NVSS (Condon et al. 1998, a catalogue of radio 
sources at 1 .4 GHz), and data from the Planck LFI and HFI. Two 
clusters in our sample exhibit relatively bright (5(1.4 GHz) > 1 
Jy), flat spectrum radio sources within a radius of 15' from the 
X-ray peak. These sources are clearly seen in LFI data and could 
potentially affect the SZ measurement. However, as we discuss 
below in Sect. 5, inclusion or exclusion of these objects has a 
negligible effect on the derived scaling relations. 



4.2.2. Impact of assumptions on pressure profile and scaling 

The blind SZ signal detection method used to detect and extract 
the ESZ clusters from the Planck survey (Planck Collaboration 
201 Id) implements the universal pressure profile from Arnaud 
et al. (2010). More specifically, the baseline model makes use of 
the generalised NFW profile fit to the 31 individual REXCESS 
cluster pressure profiles, after removal of the mass dependence 
by scaling according to the Msoo - Yx.500 relation given in Eqn. 1 . 
However, Arnaud et al. showed that the scatter of the individ- 
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ual cluster pressure profiles about the universal form increases 
toward the central regions, since cool core systems are more 
peaked, and morphologically disturbed systems are shallower, 
respectively, than the mean. In their Appendix C, Arnaud et al. 
give the best fitting GNFW model parameters for the average 
scaled profiles of the REXCESS cool core and morphologically 
disturbed subsamples. 

As our cluster sample contains both cool core and morpho- 
logically disturbed systems, it is pertinent to investigate the ef- 
fect of the baseline pressure profile assumption on the result- 
ing Fsoo values. We thus re-ran the Fsoo extraction process sep- 
arately for each object using the cool core and morphologically 
disturbed cluster profiles given in Appendix C of Arnaud et al. 
(2010). The X-ray size 6500 is kept the same in each case so 
that we are investigating the impact of the pressure profile shape 
within a fixed aperture. Figure 3 shows the ratio of the Fsoo of 
the cool core and morphologically disturbed profile extractions 
to that of the baseline model. 

Clear trends are seen in both cases: the ratio tends to increase 
(decrease) with 9sqq if the morphologically disturbed (cool core) 
profile is used instead of the baseline universal profile. Up to 
0500 ~ 10' the ratio differs from unity only by 2 per cent on aver- 
age. Beyond ^500 ~ 10', the derived Fsoo starts to differ gradually 
from the baseline value. This effect can be traced to the influence 
of the Planck angular resolution. Since the SZ signal extraction 
uses all Planck-HFl channels, the effective angular resolution 
is that of the channel with the largest FWHM (~ 10 arcmin at 
100 GHz). Below this angular scale the profile shape is washed 
out by the convolution with the Planck beam, while above it, 
clusters are increasingly well-resolved. The two panels of Fig. 3 
show that at the largest 6500 the maximum excursion is ~ 10 
per cent. Beyond 10 arcmin, the average excursions are ~ 6 and 
~ 7 per cent, respectively, for cool-core and morphologically 
disturbed profiles. Note that the effect is symmetric, in that for 
large 6500 a cool core profile and a morphologically disturbed 
profile return a value of Y500 that differs from the baseline value 
by approximately the same amount, but the former is lower and 
the latter is higher. 

In the following, the difference in Y^qq derived from extrac- 
tion with the cool core and morphologically disturbed cluster 
profiles is added in quadrature to the uncertainty on the isoo 
from the baseline extraction. We expect this conservative error 
estimate to account for any difference in the underlying pressure 
profile shape from the universal baseline model. As detailed be- 
low in Sect. 5, we have further checked the effect of the pressure 
profile assumption on the derived scaling relation fits, finding it 
to be entirely negligible. 



5. SZ scaling relations 

We fitted the parameters governing a scaling relation between 
D\ ¥500, the spherically-integrated SZ signal within Rsm, and its 
X-ray analogue Fx.soo- We also fitted parameters governing scal- 
ing relations between Fsoo and various other X-ray-derived 
quantities including Mg.500, Tx and Lx,500- We further investi- 
gated the relation between Ysqq and the total mass, M500, us- 
ing the M500 - ix,50o calibration given in Eqn. 1 . 



5.1. Fitting metliod 

For each set of observables (B,A), we fitted a power law relation 
of the form£'(z)'' Y5Q0 = 10^^ [EizY X/Xq]'^, where E(z) is the 
Hubble constant normalised to its present day value and y and 



K were fixed to their expected self-similar scalings with z. The 
fit was undertaken using linear regression in the log-log plane, 
taking the uncertainties in both variables into account, and the 
scatter was computed as described in Pratt et al. (2009). In brief, 
assuming a relation of the form Y = aX + b, and a sample of 
data points (F,, X^,) with errors ctyi and cr^, , the raw scatter was 
estimated using the error-weighted distances to the regression 
line: 



1 



N-2 



Y^Wi{Yi-aXi-bf 



where 



and 



2 22 
CTj, -I- a cr^, . 



(2) 



(3) 



The intrinsic scatter a-\ was computed from the quadratic differ- 
ence between the raw scatter and that expected from the statisti- 
cal uncertainties. 

We use the BCES regression method (Akritas & Bershady 
1996), which takes into account measurement errors in both co- 
ordinates and intrinsic scatter in the data and is widely used in 
astronomical regression, giving results that may easily be com- 
pared with other data sets fitted using the same method. We fitted 
all relations using orthogonal BCES regression. 

5.2. Effect of point sources and ctioice of baseline pressure 
profile 

We fitted the SZ scaling relations excluding the two objects with 
significant radio source contamination (Sect. 4.2. 1). For all rela- 
tions the change in normalisation, slope and intrinsic scatter is 
negligible compared to their associated uncertainties. We there- 
fore consider the contamination by radio sources to have a neg- 
ligible effect on scaling relation fits and proceed with the full 
sample of 62 clusters. 

We have also checked whether the best fitting scaling rela- 
tions are affected by the choice of baseline pressure profile, as 
follows. For the cool core subsample, we assigned the Fsoo de- 
rived from extraction using the cool core pressure profile. Of the 
remaining 40 systems, for the 20 objects with the lowest cen- 
tral density (Fig. 2), we assigned the Ksoo value derived from 
extraction using the morphologically disturbed profile. We then 
re-fitted all the scaling relations. The resulting best fits are in 
full agreement with those derived from the baseline universal 
profile; i.e., the difference in best fitting parameters (2 per cent 
maximum) is again smaller than their respective uncertainties. 

5.3. Correction for selection bias 

It is well known that scaling-relation determinations of the 
sort we are considering can be biased by selection effects of 
Malmquist and Eddington type when a significant part of the 
sample lies near a selection cut (for discussions in a cluster con- 
text see e.g., VikhHnin et al. 2009; Pratt et al. 2009; Mantz et al. 
2010 and Andersson et al. 2010). We estimated the effect of the 
Planck SZ selection as follows. In order to impose a selection 
cut on the mock catalogues, we used the observed relation be- 
tween Fgoo and S/N from the region significantly above the 
selection cut and extrapolated below it, along with an estimate 
of scatter again from observations, carried out in several red- 
shift bins. We then constructed large mock catalogues of clusters 
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Table 2: Best fitting orthogonal BCES parameters for scaling relations. 







Relation 


^ohs 




C log.i 






r 


K 






i'soo 


- c* ix.soo 


-4.02 ±0.01 


0.95 ± 0.04 


0.10 ±0.01 


-4.02 


0.96 






1 X lO"'* Mpc^ 




Ysm 


- M500 


-4.19 ±0.01 


1.74 ±0.08 


0.10 ±0.01 


-4.22 


1.74 


-2/3 




6x 10'" Mo 




isoo 


-Tx 


-4.27 ± 0.02 


2.82 ±0.18 


0.14 ±0.02 


-4.22 


2.92 


-1 




6keV 




5^500 


- Mg,500 


-4.05 ± 0.01 


1.39 ±0.06 


0.08 ± 0.01 


-4.03 


1.48 


-2/3 




1 X lO"* Mo 




5^500 


-L[0.1-2.4]x,5oo 


-4.00 ± 0.02 


1.02 ±0.07 


0.14 ±0.02 


-3.97 


1.12 


-2/3 


-112, 


7 X lO'*" erg s"' 




I'soo 


-L[0.5-2]x.5oo 


-3.79 ± 0.03 


1.02 ±0.07 


0.14 ±0.02 


-3.75 


1.12 


-2/3 


-112 


7 X lO"*"* erg s"' 



Notes, c* = (ctt / ineC')/ (figinp). 

Relations are expressed as E(zy [Dj^ Ysoo] = 10^ [E(zY X/Xq]^. The logarithmic intrinsic scatter of the relation is denoted by crio. 




0.001 0.010 10-^ 10" 10-^ 

D^' (ot/iti, c')/(He mp) Yx, 500 [arcmin'] (ot/iti^ c')/(^e nip) Yx, [Mpc'] 

Fig. 4: SZ flux vs X-ray prediction. Blue stars indicate cool core systems. Left panel: Relation plotted in units of arcmin^. The dashed line is the 
prediction from REXCESS X-ray observations (Arnaud et al. 2010). Right panel: Relation plotted in units of Mpc^. The SPT results are taken 
from Andersson et al. (2010). 



through drawing of Poisson samples from a suitably-normalised 
Jenkins et al. (2001) mass function; to each cluster we assigned 
a value of Fgoo by adopting scaling relations with scatter that 
are consistent with the observed values. This procedure leads to 
a predicted S/N value that can be used to impose selection cuts 
on the mock sample. We applied it to the full 158 cluster sample 
as the only X-ray information needed was the position for SZ 
signal re-extraction. 

The effect on scaling relations is then assessed by assigning 
further physical properties to the mock catalogue. Following the 
methods of the X-ray analysis, ix.soo is obtained directly from 
the mass using Eqn. 1, while Mg_5oo and Lx,500 are obtained from 
assumed input scaling relations including scatter. Finally Tx is 
simply obtained from Fx.soo/A^g.soo on a cluster-by-cluster basis. 
The input scaling relation slopes and amplitudes are then ad- 
justed until the mock observed samples match those recovered 
from actual data in Table 2. The input slopes then provide an es- 
timate of the bias-corrected slope that would have been obtained 
had the bias been absent. The original and bias-coiTected esti- 
mates are shown in Fig. 4 and 5, and the best-fitting parameters 
for each relation are given in Table 2. Note that the slopes of the 



Jx.soo and Tx relations are derived quantities fixed by the other 
scalings we have chosen. 



As seen in Table 2 and in Fig. 5, the importance of the 
Malmquist correction depends on the relation under considera- 
tion. In the case of Fx.soo and M500 it is negligible, due to the very 
small scatter seen in these relations. For the other relations, how- 
ever, the Malmquist corrections can be comparable to the quoted 
observational uncertainties, indicating that despite the dynamic 
range of the Planck sample, there are biases introduced by the se- 
lection cut. The bias-corrected slopes in Table 2 are thus our best 
current estimates of the true underlying slopes. One should also 
bear in mind that the bias correction itself carries uncertainty, 
which we have not been able to estimate, and this increases the 
uncertainty on the underlying slope. 



Note in particular that the bias correction leaves the F500- 
Yx,5m relation completely consistent with the expected slope of 
unity, while the relation to M500 remains consistent with a slope 
5/3. 
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Fig. 5: Scaling relations for the 62 clusters in the Planck-XMM-Newton archive sample; fits are given in Table 2. Cool core systems are plotted as 
blue stars, other systems as black dots. In the upper panels, the dotted line denotes the observed scaling relation fit. In the lower panels, the dotted 
line denotes the observed scaling relation fit, while the solid line shows the fit once the effects of selection bias are taken into account. The grey 
shaded area indicates the Icr uncertainty. The SPT results are taken from Andersson et al. (2010). 



6. Discussion 

6.1. SZ flux vs X-ray prediction 

Figure 4 shows the fundamental relation probed by the present 
study, that between the measured quantities Fx.soo and Y^qq. We 
recall that the link between these two quantities is sensitive to 
the structure in temperature and density. Note that X-ray infor- 
mation is used to determine the radius of integration for the SZ 
signal (i.e., R^qq) and its overall shape (i.e., the underlying uni- 
versal pressure profile). However, as we have shown above in 
Sect. 4.2.2, the amplitude of the SZ signal is relatively insen- 
sitive to the assumed pressure profile shape, so that the use of 



X-ray priors reduces to a choice of integration aperture. Thus 
we regard the X-ray and SZ quantities as quasi-independent. 

In the left-hand panel of Fig. 4 the relation is plotted in units 
of arcmin^, and shows the excellent agreement between the ob- 
served Igoo, lx,500 and the X-ray prediction from REXCESS 
(dashed line). Indeed, fitting the relation with the slope fixed to 
unity yields a normalisation Fgoo/ Jx,500 = 0.95 + 0.03, perfectly 
consistent with the value 0.924 ± 0.004 found for the REXCESS 
sample (Arnaud et al. 2010), and less than unity as expected for 
radially-decreasing temperature profiles (e.g., Pratt et al. 2007). 
Furthermore, the relation is quite tight (see below), and there is 
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no indication that cool core systems differ systematically from 
the other systems. 

The right-hand panel of Fig. 4 shows the relation be- 
tween Fx,50o and the spherically-integrated Compton parameter 
D\ F500 ■ Note in particular that the slope of the bias-corrected re- 
lation is completely consistent with unity, and that the intrinsic 
scatter (calculated as described above in Eqns. 2 and 3) is excep- 
tionally small, cTjog i = 0.09 + 0.01. A comparison with recent 
results obtained by SPT shows a slight difference in normalisa- 
tion, although it is not significant given the larger uncertainties 
in the latter measurement (Andersson et al. 2010). 

6.2. Scaling relations 

In this Section we investigate other relations between F500 
and X-ray quantities. Note that since M500 is derived from ix.soo, 
its dependence on Fsoo is directly linked to the F500 - 
Yx,5Q0 relation with the exception of differing E(z) dependencies. 
Moreover, Mg ^oo and Tx are not independent; they are related 
via Eqn. 1 . It is still useful to investigate these relations, though, 
both for completeness and for comparison to recent results from 
ground-based studies. 

Relations between Fgoo and gas mass Mc, 500 and the X- 
ray temperature Tx are shown in the upper panels of Fig. 5. The 
bias-corrected relations yield slopes that are consistent with self- 
similar (5/3 and 5/2, respectively) to high accuracy. Scatter about 
the Fsoo - Mg.soo relation is small, at criog j - 0.09 + 0.01, 
while that about the Fsoo - Tx relation is among the largest 
of the relations, at criog j = 0.14 + 0.02. Once again, cool core 
systems are fully consistent with the other clusters, and there is 
no particular evidence that cool core systems show less scatter 
than the sample as a whole. 

The bottom left-hand panel of Fig. 5 shows the relation 
between isoo and mass. Here again, the slope of the bias- 
corrected relation is fully consistent with self-similar (5/3). The 
scatter is small (criog j - 0.10 + 0.01), although it is a lower limit 
since the scatter between lx.500 and total mass is not taken into 
account in deriving M500 (it is in fact the same as that about 
the F500 - Jx,500 relation except for the different E(z) scal- 
ing). Both slope and normalisation are in excellent agreement 
with X-ray predictions from REXCESS, as expected from the 
good agreement in the Y^qq - Fx.soo relation. The slight off- 
set in normalisation of the relation found by SPT (Andersson 
et al. 2010) can be explained by the different calibration of the 
Msoo - Yx,50o relation used in their study; it is not a significant 
offset given their larger normalistion uncertainties. 

The ease of detecting clusters through their X-ray emission 
makes the X-ray luminosity an important quantity, and its cali- 
bration with the SZ signal is imperative for maximising the syn- 
ergy between the Planck all-sky survey and previous all-sky X- 
ray surveys such as the RASS and the upcoming eROSITA sur- 
vey. The slope of the D\ Y^qq - Lx.soo relation for the present 
sample, 1.12 + 0.08, is in excellent agreement with the slope pre- 
dicted from X-ray observations alone (1.07 + 0.08, Arnaud et al. 
2010 from REXCESS), and the normalisation is also consistent 
within the uncertainties. The slight offset in the best fitting nor- 
malisation for the present sample relative to the REXCESS pre- 
diction can be attributed to the relative lack of strong cooling 
core clusters in the present sample compared to REXCESS (see 
Fig. 2). The scatter, cr\og i - 0.14 + 0.02, is largest about this 
relation due to the influence of cool cores, which are segregated 
from the other systems and all lie to the high-luminosity side. 
Indeed, as Fig. 5 shows and Table 2 quantifies, while the vast 



majority of the dispersion about the F500 - ix,500 relation is 
due to cool cores, these systems do not contribute significantly to 
the dispersion about the Y^qq - Fx.soo relation. Thus while the 
X-ray luminosity is very sensitive to the presence of cool cores, 
D\ Fgoo appears to be less so. 

The slope of our best fitting Fgoo - Lx.soo relation is 
also fully consistent within Icr with that derived by Planck 
Collaboration (201 If), which is based on a bin-averaging anal- 
ysis at the position of known X-ray clusters in the MCXC 
(Piffaretti et al. 2010). As X-ray selection is more sensitive to 
the presence of cool cores (due to the density squared depen- 
dence of the X-ray luminosity), one might expect the Planck 
Collaboration (201 If) best-fitting Dj^Y^oo - Lx^soo relation to 
be shifted to slightly higher luminosities (i.e., a slightly lower 
normalisation), as is seen. However, since the MCXC selection 
function is both complex and unknown, it is impossible to cor- 
rect their relation for Malmquist bias effects. Thus some part of 
the normalisation difference between the two studies arises from 
correction for selection effects. However the good agreement be- 
tween our results and those from Planck Collaboration (201 lf)'s 
bin-averaged analysis argues that sample selection does not have 
a strong effect on the scaling relations derived from the latter 
analysis. 

Finally, the results we have derived show that X-ray and 
SZ measurements give a fully coherent view of cluster struc- 
ture out to moderately large scales. In particular, they indicate 
that effects of clumping in the X-ray gas are not significant, at 
least in the mass and radial range we have probed in the present 
study. Furthermore, the excellent agreement between the ob- 
served Fjoo - lx.500 relation and the X-ray predictions argue 
that the SZ and X-ray calibrations we have used are fundamen- 
tally sound. 



7. Conclusions 

We have presented SZ and X-ray data from a sample of 62 lo- 
cal (z < 0.5) galaxy clusters detected at high S/N in the Planck 
survey and observed by XMM-Newton. The objects range over 
approximately a decade in mass (M500 ~ 2 -20 x lO'"* M©), and, 
while the sample is neither representative nor complete, it con- 
stitutes the largest, highest-quality SZ-X-ray data set currently 
available. This study has been undertaken in the framework of a 
series of papers on cluster astrophysics based on the first Planck 
all-sky data set (Planck Collaboration 201 ld,e,f,h). 

SZ and X-ray quantities have been extracted within Rsoq and 
we have presented a detailed study of the resulting SZ scaling 
relations. Moreover, we have investigated how selection effects 
influence the results of the scaling relation fits. Their influence 
is subtle, but the slopes and normalisations of the scaling re- 
lations are generally in good agreement with X-ray predictions 
and other results after accounting for the selection effects. For 
the fundamental Fgoo - Ix.soo relation, we measure a remark- 
ably small logarithmic intrinsic scatter of only (10 + 1) per cent, 
consistent with the idea that both quantities are low-scatter mass 
proxies. 

The results are fully consistent with the predictions from X- 
ray observations (Arnaud et al. 2010) and with recent measure- 
ments from a smaller sample spanning a wider redshift range 
observed with SPT (Andersson et al. 2010). The results are also 
in excellent agreement with the statistical analysis undertaken 
at the positions of known X-ray clusters (Planck Collaboration 
201 If). This excellent agreement between observed SZ quanti- 
ties and X-ray-based predictions underlines the robustness and 
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consistency of our overall view of ICM properties. It is diffi- 
cult to reconcile with the claim, based on a recent WMAP 7-year 
analysis, that X-ray data over-predict the SZ signal (Komatsu 
et al. 2011). 

The results presented here, derived from only 62 systems, 
provide a maximally-robust local reference for evolution stud- 
ies or for the use of SZ clusters for cosmology. Overall, the 
agreement between the present results, ground-based results and 
X-ray predictions augurs well for our understanding of cluster 
astrophysics and for the use of clusters for cosmology. Future 
work will make use of the individual pressure profile shape as 
derived from X-rays to further improve the SZ flux extraction. 
Comparison of X-ray and SZ pressure profiles will also be un- 
dertaken, as will comparison of measurements with independent 
mass estimation methods. We will also extend our analysis to the 
full Planck catalogue, observing higher-redshift systems, to con- 
strain evolution, and lower-mass objects, to better probe cluster 
astrophysics. 
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Appendix A: Optimised SZ extraction and 
comparison with X-ray predictions 

As discussed in the main text, with the present cluster sample we 
have optimised the SZ photometry by using the higher-quality 
estimate of the X-ray size 6500, derived from /?goo measured us- 
ing the Msoo - Yx,50Q relation as detailed in Sect. 3.3. 

In Fig. A. 1 we examine the change in F500 when derived us- 
ing different characteristic sizes 6500 to extract the SZ signal. We 
also compare the SZ signal predicted using X-ray observations 
(from the Fsoo/^x.soo relation of Arnaud et al. 2010) to the ob- 
served SZ signal. In all cases, the ICM pressure is assumed to 
follow the baseline universal profile of Arnaud et al. (2010). As 
extensively described in Planck Collaboration (201 Id), the SZ 
flux Fsoo is computed by integrating along the line-of-sight and 
normalising the universal pressure profile. Each profile is trun- 
cated at 5 X R500, effectively giving a measure of the flux within 
a cylinder of aperture radius 5 x R500, and then converted to the 
value in a sphere of radius Rsoo for direct comparison with the 
X-ray prediction. 

The left hand panel (a) of Fig. A.l shows the relation be- 
tween Fsoo and that predicted from XMM-Newton observations 
as used in the present paper, illustrating the tight agreement be- 
tween the two quantities. We recall that here, the /?5oo within 
which the SZ signal is extracted is derived from the measured 
5'x,500 using the M500 - ix,500 relation given in Eqn. 1 . 

The middle panel (b) of Fig. A. 1 assumes that only the X-ray 
position and luminosity of the cluster are known. In this case the 
mass is derived from the Mgoo - Lx,500 relation as described in 
the MCXC of Piffaretti et al. (2010), thus yielding the charac- 
teristic size used to extract the SZ signal, 6500. This mass is also 
used to predict Fx,500 via the M500 - J'x,500 relation in Eqn. 1. 
Consistently, the expected SZ signal is extracted from a region 
of size 0500 centred on the X-ray position given in MCXC (as in 
Planck Coflaboration 201 If). 

In the right hand panel (c) of Fig. A. 1 the position and size 
of the cluster are unknown, thus Fgoo is devired blindly together 
with the SZ flux. However, the predicted SZ flux is derived as 
above in panel (b). The agreement between measured and pre- 
dicted values clearly degrades dramatically from panels (a) to 
(c). Comparing panels (b) and (a), there is a systematic shift to 



lower predicted }x,500 values, with a segregation now appear- 
ing between cool cores and the other systems. This can be ex- 
plained by the fact that using the luminosity as a simple mass 
proxy leads to an underestimate of the mass for morphologically 
disturbed systems in view of their position with respect to the 
mean Lx,500 - Msm relation (Pratt et al. 2009). The inverse effect 
is seen for the cool cores. In addition, there is a smaller impact 
on the measured Ysm via the effect of the assumed 6500 • However 
the effect is smaller: the average ratio of XMM-Newton and 
MCXC characteristic sizes 05oo,mcxc/S5OO,xmm is 0.95 + 0.06, 
corresponding to a change in area of ~ 10 per cent, which trans- 
lates into a similar variation in SZ flux. This shows that the X-ray 
luminosity in the MCXC is a sufficiently good mass proxy for a 
reliable size estimate. 

Finally, panel (c) of Fig. A. 1 illustrates the size-flux degen- 
eracy in blind Planck measurements. When Fsoo is measured 
blindly, the size is on average overestimated (see also Planck 
Collaboration 201 Id), and so the disagreement with predictions 
is even more apparent. 

We see that as a result of the size-flux degeneracy, an ac- 
curate estimate of the characteristic size is mandatory in or- 
der to derive an accurate measure of Fgoo- A similar conclusion 
was reached in Planck Collaboration (201 Id, see their Fig. 11), 
where the effect was demonstrated using the full sample of 158 
clusters known in X-rays (i.e., those included in the MCXC). 
These authors found that, in addition to a reduction in intrinsic 
scatter (from 43 to 34 per cent), knowledge of the cluster size 
dramatically reduced the offset of the measured Y^qq to that pre- 
dicted from X-rays (from 84 to 14 per cent - compare panels (c) 
and (b) of Fig. A. 1 above). 

However as noted in Planck Collaboration (201 Id), there is 
still a small but systematic discrepancy. This is mostly due to the 
use of Lx,50o, a quantity which shows considerable scatter with 
mass, as a mass proxy. The superior constraints provided by the 
XMM-Newton observations on the cluster size and on }x,50o sup- 
press most of this remaining systematic effect (compare panels 
(b) and (a)). A smaller contribution is liked to effects due to the 
nature of the sample selection. This illustrates that a fully coher- 
ent approach is needed when undertaking a proper comparison 
between SZ and X-ray predictions. 
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Fig. A. 1: Comparison of the Wane/: measured SZ flux, Ysoo and the predictions from X-ray measurements, (crT/insC^)/(jicntp)D^Yx.5Qo- Cool core 
systems are marked as blue stars, other systems as black dots. (Left) Kx.soo from XMM-Newton (see top-left panel of Fig. 5), and Ijoo computed at 
the position and 6500 derived from XMM-Newton measurements. (Middle) Yx.500 and Ksoo respectively from Lx.500, and position and 6500 as given 
in the MCXC (Piffaretti et al. 2010). (Right) Kx.soo from Lx.soo as given in the MCXC, while Y^oo are blind Planck measurements. 
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